Skyrmion quantum numbers and quantized pumping in two dimensional topological chiral magnets 
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We investigate the general conditions to achieve the adiabatic charge and spin polarizations and quantized 
pumping in 2D magnetic insulators possessing inhomogeneous spin structures. In particular, we focus on the 
chiral ferrimagnetic insulators which are generated via spontaneous symmetry breaking from correlated two 
dimensional topological insulators. Adiabatic deformation of the inhomogeneous spin structure generates the 
spin gauge flux, which induces adiabatic charge and spin polarization currents. The unit pumped charge/spin 
are determined by the product of two topological invariants which are defined in momentum and real spaces, 
respectively. The same topological invariants determine the charge and spin quantum numbers of skyrmion 
textures. It is found that in noncentrosymmetric systems, a new topological phase, dubbed the topological chiral 
magnetic insulator, exists in which a skyrmion defect is a spin- 1/2 fermion with electric charge e. Considering 
the adiabatic current responses of generic inhomogeneous systems, it is shown that the quantized topological 
response of chiral magnetic insulators is endowed with the second Chem number. 
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I. INTRODUCTION 

After the seminal papers written by Jackiw and Rebbi' and 
Su, Schrieffer, and Heeger^, quantum number fractionaliza- 
tion has become one of the most fundamental and fascinating 
concepts in condensed matter physics. For one dimensional 
spinless fermions, a domain wall soliton induces a localized 
mid-gap state which has a half of the electric charge. Later 
this idea is further generalized by Goldstone and Wilczek.^ 
Considering slow variations of the background mass fields in 
space and time, adiabatically induced currents are derived. 
The time component of the adiabatic current describes the 
charge density induced by the spatial variation of the mass 
soliton. The spatial integration of this charge density gives rise 
to the half electric charge.^ Moreover, the temporal variation 
of the background mass fields generates the spatial component 
of the adiabatic current, which describes the adiabatic charge 
polarization or even the quantized charge pumping.** There- 
fore the existence of the topological soliton with nontrivial 
quantum numbers reflects the quantized dynamical response 
of the system. 

The discovery of the time-reversal invariant spin Hall insu- 
lators'' has triggered great attentions to the two dimensional 
topological band insulators (TBIs). TBIs are the insulating 
phases that carry gapless edge states at the sample boundary 
due to the nontrivial bulk Berry phase, which can either break^ 
or preserve the time reversal symmetry.^ Because of the dis- 
tinct bulk topological properties, tt flux defects in TBIs carry 
nonzero spin or charge quantum numbers leading to Jackiw- 
Rebbi type topological solitons in two dimensional systemsi^ 
In addition, in strongly interacting electron systems, it is pro- 
posed that TBIs can be obtained from the spontaneous break- 
ing of spin rotation symmetryi^ Since this topological "Mott" 
insulator possesses a vector order parameter N, skyrmion- 
type topological solitons are allowed, which are also expected 
to have nontrivial quantum numbers. For example, in the case 
of the spin Hall "Mott" insulator, a skyrmion defect is a boson 



carrying a quantized electric charge 2ei^ 

In this work, we focus on a different class of topologi- 
cal Mott insulators, the chiral magnetic insulators (CMIs), 
which can be obtained by considering electron-electron inter- 
actions in TBIs. In CMIs, the vector order parameter N is 
given by the stagggered spin ordering, which can acquire ex- 
plicit time dependence by coupling to the external magnetic 
field. Because of the quantized charge and spin Hall effects, a 
skyrmion defect in CMIs can carry charge and spin quantum 
numbers satisfying fermionic or bosonic statistics depending 
on the magnitude of the bulk topological invariants. Consider- 
ing smooth variations of the spin ordering directions in space 
and time, reminiscent of the approach pursued by Goldstone 
and Wilczek, we have shown that the quantized charge and 
spin pumpings are possible under the adiabatic change of the 
inhomogeneous spin structures such as the domain wall with 
rotating spins in collinear ferrimagnets or the conical ferri- 
magnetic phases in which the spatial modulation of the spin 
direction spreads over the whole system. 

In CMIs, the quantum of the pumped charge/spin or the 
quantum number of a skyrmion strongly depends on the na- 
ture of the parent TBIs. In the case of the CMIs generated 
from the anomalous Hall insulator (AHI), the quantum of the 
pumped charge/spin is solely determined by the spin/charge 
Chern numbers of the collinear magnetic ground states. When 
the magnitude of the spin order parameter is small, the charge 
and spin Chern numbers are the same as those of AHI. In 
this case, the spin can be pumped in the unit of h while 
the charge pumping is forbidden, which, at the same time, 
means that the skyrmion is a chargeless spin-1 boson. As 
the magnitude of the ordered spin moment increases, a new 
CMI phase dubbed the topological chiral magnetic insulator 
(TCMI) emerges in noncentrosymmetric systems, which sup- 
ports a fermionic skyrmion with charge e and spin |. Here 
the quantized charge and spin pumpings can be achieved. 

On the other hand, in the case of CMIs derived from the 
spin Hall insulator (SHI), the spin anisotropy inherent to SHI 
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due to the spin-orbit coupling imposes constraints on the adi- 
abatic pumping and skyrmion quantum numbers. The polar- 
izations strongly depend on the relative orientation of the spin 
anisotropy and the magnetic ordering directions. When the 
spin order parameter is parallel to the spin anisotropy direc- 
tion, the charge polarization induced by the rotating spins in 
domain walls is always zero. On the other hand, if the spin or- 
dering is perpendicular to the spin anisotropy direction, finite 
charge polarization is expected for the Neel domain walls and 
the transverse conical spin states. In the case of a skyrmion, 
it does not carry a charge quantum number independent of the 
magnitude of the spin order parameter. However, a meron de- 
fect can be charged but the net charge is not quantized. 

In addition, it is shown that the charge and spin polariza- 
tions in CMIs can be understood in the framework of the 
adiabatic topological responses of generic inhomogeneous 
crystalsiii^ The inhomogeneity induced polarization current 
contains a topological part whose topological nature is en- 
dowed with the second Chern number, which is quantized in 
four dimensional closed manifolds. When the inhomogeneity 
is introduced by the spatially varying three-component unit 
vector N, the inhomogeneity induced topological current can 
be separated into the homogeneous and inhomogeneous parts, 
which allows the second Chern number to be represented by 
the product of two independent topological invariants. One is 
the Chern number that is quantized in the momentum space 
and the other is the skyrmion number defined in the posi- 
tion/time spaces. 

Since electron correlation effects are required to achieve 
CMI in addition to large spin-orbit coupling, transition metal 
oxides with 5d electrons can be an ideal platform to re- 
alize CMI. Recently, it is confirmed that the ground state 
of Sr2lr04 is a spin-orbit entangled Mott insulator with 
anti-ferromagnetic spin ordering,i^— which manifests the 
strong interplay of electron correlation and spin-orbit cou- 
pling effects in Ir-based compounds. In the case of Na2lr03 
which is expected to be SHI in non-interacting limit, anti- 
ferromagnetic spin ordering develops when electron correla- 
tion effects are included on the parent SHI phase.'** In this 
respect, Na2lr03 is a candidate material in which the topolog- 
ical response of CMI can be observed. Moreover, it is recently 
proposed that bilayers of the perovskite-type transition metal 
oxides grown along the [111] direction can be potential can- 
didates for two-dimensional topological insulators.''''^" The 
transition metal ions in the bilayer are located on a honeycomb 
lattice consisting of two trigonal sublattices on different lay- 
ers. Interestingly, since the layer potential difference can be 
easily created by applying an electric field or by sandwiching 
the bilayer between two different substrates, lattice inversion 
symmetry can be explicitly broken in this system. Therefore 
to reveal the ground state of the bilayer systems, it is important 
to understand the interplay of the spin-orbit coupling, electron 
correlation, and the broken inversion symmetry. 

The rest of the paper is organized in the following way. In 
Sec. |II]we construct a model Hamiltonian of CMIs taking into 
account electron-electron interactions in TBIs on the honey- 
comb lattice. Sec. |III]we consider the general conditions to 
achieve the quantized charge and spin pumpings through a 



domain wall for CMI derived from AHI. The charge and spin 
quantum numbers of a skyrmion in CMI are also discussed. 
The investigation of the inhomogeneity induced charge polar- 
ization and skyrmion/meron quantum numbers is extended to 
CMI dervied from SHI in Sec.|lV] In Sec.[V]the adiabatic po- 
larization currents for generic inhomogeneous systems are de- 
rived by using the semi-classical gradient expansion method. 
The topological origin of the quantized charge pumping in 
CMI is discussed from this viewpoint. Finally, we conclude in 
Sec. [Vl] The details of the procedures to derive the adiabatic 
polarization currents from the gradient expansion method are 
presented in Appendix lAl 

II. LATTICE MODEL OF TOPOLOGICAL CHIRAL 
MAGNETIC INSULATORS 

The magnetic insulator with broken inversion symmetry 
is an ideal playground in which inhomogeneity induced re- 
sponses can be investigated. This is because the lack of the 
spatial inversion symmetry allows Dzyaloshinskii-Moriya in- 
teraction which leads to spatially modulated spin structures. 
In this work, we focus on the magnetic insulators which 
can be obtained spontaneously from TBIs including inver- 
sion symmetry breaking terms. To demonstrate the main idea 
concretely, we start by constructing a model Hamiltonian of 
TCMI on the honeycomb lattice. Up to now, two types of TBIs 
are proposed on the honeycomb lattice. One is the anoma- 
lous Hall insulator (AHI) with broken time-reversal symmetry 
and the other is the time-reversal invariant spin Hall insulator 
(SHI). Starting from these TBIs, magnetic ordering can be de- 
veloped when electron correlation is considered;^ 

The Hamiltonian describing the TCMI can be written in the 
following way, 

H = Hi + H2 + Hu, (1) 

in which 

where (ij) and {{ij)) denote the nearest-neighbor (NN) and 
next-nearest-neighbor (NNN) pairs, respectively, ti indicates 
the hopping amplitude between NN sites. t2 (cr) is a spin inde- 
pendent constant t2{(j) = t2 in AHI while ^2(0') = ^2 (— ^2) 
for spin-up (spin-down) electrons in SHI, where t2 is the am- 
plitude of the pure imaginary hopping between NNN sites. 
The unit vector ey is defined as etj — {dlj x d^^ ) / 1 d^^ ^ ^Ijl 
where dlj and dfj are the bond unit vectors along the two 
bonds which are traversed by the electron when it moves from 
the site j to i. The noninteracting Hamiltonian Hi has a semi- 
metallic ground state at half-filling with two Dirac points at 
the momenta Ki and K2(= — Ki). The staggered chemical 
potential ^s.i in H2, which breaks the lattice inversion sym- 
metry, satisfies ^s,i = Ms (~Ms) if i belongs to the A (B) 
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FIG. 1 : (Color online) The schematic phase diagram of the full lattice 
mean field Hamiltonian -&fu11 (Eq. ((S}) for coUinear magnetic phases 
derived from AHI. (a) The phase diagram shown in two dimensional 
space of (mAF, ^^s) for given t2 > 0. (b) One dimensional section of 
the 2D phase diagram, which is obtained by varying ttiaf for given 
i2 > > along the long arrow in (a). In both (a) and (b), a dotted 
line indicates the gapless phase boundary between gapped phases. 
In (b), each phase is characterized by the charge (C'c) and the spin 
(Cg) Chern numbers. Here TCMI (SDW) means a topological chiral 
magnetic insulator (spin density wave insulator). 



sublattice. Considering the potential realization of CMI in 
bilayers of the perovskite-type transition metal oxides grown 
along the [111] direction, the inclusion of the staggered chem- 
ical potential iis,i is required. Interestingly, it is shown 
later that the existence of the inversion symmetry breaking 
term ^sa plays an essential role to stabilize TCMIs. 

H\j indicates electron-electron interactions generating a 
magnetic ground state. We introduce a collinear spin ordering 
along direction represented by Sa ~ and Sb ~ msf 
at the two sites A and _B in a unit cell. Assuming the transla- 
tional invariance of the magnetic ground state, the mean field 
approximation for H\] leads to the following Hamiltonian, 

Hm¥ = X! X! ^^(k)[mF,,;(ToCri) + mAF,i('^zO'i)]V'(k), 
k i—x,y,z 

= ^ V^^(k)[(mF • a) + r,(mAF • o')]V'(k), (2) 

k 

where V'^(k) = (c^^^(k), _^(k), 4^^(k), ^^^(k)). mF = 
rrip-r and ihaf = mApf with tof — {rriA + ■mB)/2 and 
mAF = [iTT-A — ™b)/2. CTi and (i ~ x,y,z) are Pauli 
matrices for spin and sublattice degrees of freedom and (Tq, 
To are identity matrices. Here <TiTj = ai (S) Tj indicates a 
4x4 matrix acting on the 4 component vector ^/'(k). Also 
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FIG. 2: (Color online) Evolution of the band structure as rriAF 
changes for t2 = 2, fis ~ 1. The AHI (rriAF ~ 0.5) turns into the 
TCMI (mAF = 1.5) via a critical point at ttiaf = ^2 — /is = 1. Due 
to the broken time reversal and inversion symmetries, a band touch- 
ing occurs only at k = K2 between spin-up bands, which changes 
the charge and spin Chern numbers simultaneously. 



when the direct product of Pauli matrices is involved with an 
identity matrix such as tq or (Tq, we omit the identity matrix 
in the product if the dimensionality of the Hamiltonian is ob- 
vious from the context. For example, Toai = (Ti in Eq. (|2|. 
Assuming |mAF| ^ ItjfI, we first neglect the ferromagnetic 
component tof, whose influence on the adiabatic polarization 
is considered later. 

The full lattice mean field Hamiltonian can be written as 



Hfuu — Hi + H2 + Hi 



MF, 



(3) 



which leads to the phase diagram shown in Fig. [T] and Fig. [3] 
The phase diagram is obtained in the following way. Since the 
(topological) phase transition between neighboring insulators 
occurs via band touching at k = Ki or K2, the phase diagram 
can be investigated by constructing the effective Hamiltonain 
describing the low energy particles near the two nodes k = Ki 
or K2. Explicitly, it is given by 



-fflow(q) = qxTxl^z + qyTy 



(4) 



where cr^ 
trices i'. 



x,y,z 



: do (cr^ = crj for AHI (SHI). The Pauli ma- 
are introduced to indicate the node degrees of 
freedom, i.e., Ki — 1) and K2 (i^z = ^ !)■ The mo- 
mentum {qx, Qy) is measured with respect to the nodes. When 
the Hamiltonian can be written as a massive Dirac Hamilto- 
nian such as iJ = q.j.,Tx + qyTy + rriTz, the sign reversal of the 
mass m changes the Chern number of the occupied band by 
AC = sgn(m) via a band touching at to = 0. Therefore the 
location of the phase boundary and the consequential change 
of Chern numbers can be understood by comparing the rela- 
tive magnitude of /^s, ^2, and toaf in Eq. (01). In addition, if 
the Chern number of one phase is known, the Chern numbers 
of all the other phases can be determined by investigating the 
sign change of the mass term of the effective Dirac Hamilto- 
nian. For example, if toaf ^ ^2 and toaf ^ /^s, the corre- 
sponding gapped phase, i.e., the SDW phase, should have zero 
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Chem number since it is a topologically trivial. In this way, 
the phase diagrams in Fig.[T]and Fig.|3]can be obtained. 

In addition, to confirm the structure of the phase diagram 
and topological property of each phase, the Chern number of 
each phase is numerically computed including the full band 
dispersion in the Brillouin zone. The Chern number Ca of the 
occupied band with spin a is defined as. 



(5) 



where the momentum space Berry curvature ^^(k) is de- 
fined as -Fcr(k) = dk^Ay^^{k) - dk^A.^^^{k) in which 
the Berry potential A^.o-(k) is given by ^^ ^(k) — 
-i($<,(k)|(9fc^|$<,(k)).2i'22 Here |$^(k)) indicates the peri- 
odic part of the Bloch wave function corresponding to the oc- 
cupied state with spin a. 

Interestingly, various insulating magnetic phases with dis- 
tinct topological properties are possible depending on the rel- 
ative magnitudes of ttiaf, t2, and /is. Different insulating 
phases are distinguished based on the Chern numbers obtained 
for the collinear magnetic ground states. Fig. [T] shows the 
magnetic ground states derived from AHI. Here we continue 
to use the term "AHI" as long as the magnetic phase possesses 
the charge Chern number Cc = ±2 and the spin Chern num- 
ber Cs — 0. The charge (Cc) and spin (C5) Chern numbers 
are defined as Cc — C^ + Ci and Cs — — C^. In particu- 
lar, for ||i2| - 1^5 1 1 < ^afI < ||t2| + |M5||,theTCMI,which 
can be characterized by the odd integer charge and spin Chern 
numbers ( I Cc I = \Cs\ — 1), is obtained. SDW indicates the 
magnetic insulator with Cc ~ Cs ~ 0. 

The band structures of two topological phases, i.e., AHI and 
TCMI, and the band touching at the phase boundary between 
them are shown in Fig. |2] It is interesting to note that due to 
the broken time reversal and inversion symmetries, the band 
touching occurs only between spin-up bands at k = K2. Since 
the Chern number for the occupied spin-up band changes by 
1 (AC^ = 1), the charge and spin Chern numbers also change 
simultaneously leading to TCMI phase. 

In the case of magnetic insulators derived from SHI, the 
nature of the ground state depends crucially on the relative 
orientation between mAp and t2 = t2Z. If mAF||t2, i.e., 
mAF = niApz, the phase diagram is shown in Fig. [3] (a) and 
(c). In this case, similar to the case of magnetic insulators de- 
rived from AHI, various topological insulators can be obtained 
depending on the relative magnitude of t2, /is, and toaf- Here 
"SHI" also includes the magnetic phases which have the same 
Chern numbers as the usual nonmagnetic SHI phase assuming 
that the spin z component Sz is conserved. However, if Sz is 
not conserved, it is equivalent to a topologically trivial phase. 
As the relative angle between mAF and t2 increases, the phase 
boundaries between different gapped phases change. In par- 
ticular, when mAF-Lt2, the area of TCMI phase shrinks to zero 
and the phase diagram contains only two phases, i.e., SHI and 
SDW. This is because when mAF-Lt2, in contrast to the case of 
niAF||t2, gap-closing appears at the two nodes simultaneously 
along the phase boundary between SHI and SDW. In this case, 
the Chern number of the system does not change. The varia- 
tion of the phase diagram as the relative angle between mAF 
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FIG. 3: (Color online) The schematic phase diagram of the full lattice 
mean field Hamiltonian f/puii (Eq. ([3}) for collinear magnetic phases 
derived from SHI. (a) The phase diagram shown in two dimensional 
space of (ttiaf, /is) for given t2 > when t2 || mAp. (b) The phase 
diagram for t2 -L mAF. (c) One dimensional section of the 2D phase 
diagram, which is obtained by varying rriAF for given t2 > fis > 
along the long arrow in (a). In (a), (b), (c), a dotted line indicates the 
gapless phase boundary between two gapped phases. 



and t2 changes implies that starting from a TCMI phase with 
given t2, /is, and mAF, if we rotate the direction of mAF rela- 
tive to t2, topological phase transitions should occur, through 
which TCMI turns into either SDW or SHI phase. 

From the collinear magnetic states, the inhomogeneous spin 
structures can be introduced by considering smooth deforma- 
tion of spin ordering directions. When the inversion sym- 
metry breaking term fis is much smaller than the spin ex- 
change couplings, the spins are aligned almost collinearly. In 
this limit, the local inhomogeneous spin structures, such as 
skyrmion defects and domain walls, can be developed. On 
the other hand, since the Dzyaloshinskii-Moriya interaction is 
allowed between localized spins when the inversion symme- 
try breaking term is finite, the inhomogeneous spin structure 
can be developed globally leading to spiral (or conical) spin 
orderings. All these inhomogeneous spin structures can be de- 
scribed by replacing mApf ■ a with mAF<^ ■ N{r, rf) where the 
three-component unit vector N is given by 

7V(r,r/) = (N^,Ny,Nz) 

= {svL\9{rj) cos 0(r), s\w9{rf) sin (/)(r), cos 0{ri)). 

(6) 

For example, for mAF = rn^pz, a static configuration with 
d{ri) = 9{\r\) satisfying the boundary conditions 0(|r I = 0) = 
TT and 6'(|r| — > 00) — 0, describes a skyrmion. Adopting 
if){r) — tan^^(|-), a skyrmion can be specified by the two pa- 
rameters (|r|, (p) which is nothing but the polar coordinates for 
two dimensional space. On the other hand, taking (f){r) = Q-r 
with proper boundary conditions, N{r, -q) describes a spiral 
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FIG. 4: (Color online) An example of the adiabatic variation of the 
inhomogeneous spin structure. A spiral spin state with the modula- 
tion wave vector Q = is described. A pair of the spins at the two 
sublattice sites A and B are aligned in the opposite directions. The 
magnitude of the uniform component of the rotating spins is varying 
adiabatically (along the vertical direction in the figure), which can be 
controlled in experiment by applying external magnetic fields. 



spin ordering with the modulation wave vector Q or a smooth 
variation of spin directions around a domain wall.^^ In these 
cases, the spatial modulation itself can be described by a sin- 
gle parameter. For example, taking Q — Qy, (f){r) — (t){y) 
depends only on the y-coordinate. If we adopt the time coor- 
dinate t as f], N{y, t) describes the temporal variation of the 
inhomogeneous spin structures. Then = cos 9{t) denotes 
the uniform spin component independent of the spatial coor- 
dinates and the adiabatic variation of 9{t) controls the relative 
magnitude of the rotating and uniform spin components. In 
Fig. m we describe the adiabatic temporal change of the spi- 
ral ferrimagnetic ordering. Here the magnitude of the uniform 
component decreases as a function of t. It is the main purpose 
of this paper to understand the topological responses of the 
system due to the adiabatic change of the N which varies in 
the two dimensional parameter spaces, either (|r|, 0) or (y, t). 



III. TOPOLOGICAL RESPONSES OF CHIRAL 
MAGNETIC INSULATORS DERIVED FROM AHI 

A. Effective action for the charge and spin responses 

Let us first consider the charge and spin responses of CMIs 
which are derived from AHI. The low energy Hamiltonian, 
which is obtained by linearizing the energy spectrum of ifpuii 
near the two Dirac nodes, can be written as 

i/eff = J (fri)''{r,t)[{~ida:)T^Vz + {-idy)Ty 

+ MsT^ + t2T:,Vz -h mAFrziV(r, t) ■ a]ilj{r, t), (7) 



where the Fermi velocity of the Dirac particle is scaled to 1 
and h is set to 1 . The corresponding action is given by 

j dtd^rij^{r,t)-f'[{-ird^,) 

- ifJ-s - it2Vz ~ imApNir, t) ■ CT]?A(r, t), (8) 

where = t,x,y) is defined as (7*, 7^, 7^) — 

{—iTz, TyVz, —Tx) and the summation over repeated indices 
is assumed throughout the paper. 

We consider a local unitary transformation U satisfying 
U^{a ■ N)U = Uz, which rotates each spin to the +z direc- 
tion. Defining il) = Uip', the effective action can be written 
as 

S=J dtd'r^'^ir,t)^'[ri-zd, + A';:-^a-t^) 

- ifis ~ ihi^z - imAF(TzWir, t), (9) 

where ■ il^ = iU\d^U), which shows that the inhomo- 
geneous spin structure induces SU(2) spin gauge fields in the 
rotated frame. The electromagnetic U(l) gauge fields A'^ are 
introduced via minimal coupling and the electron charge — e 
is set to -1. To derive the effective action for gauge fields, we 
split the action into two pieces S' = S*! + S'2 in which 

/d^k 

f d^k d^q , 

X A^(a,^Go i)M^(g) - • t ,{q)]i^' (k) , 

(10) 

where 

G(7^(fc) = CJ - [kxT^Vz + kyTy + ^sTz + t2Tzl'z + mAfTzCTz] 

(11) 

Here the three momenta k is defined as k — {uj,kx,ky) and 
(AtjAjjjAy) = (—1,1,1). The effective action can be ob- 
tained after integrating out the fermion fields and expanding 
the resulting action in powers of the gauge fields and their gra- 
dients. Straightforward calculation gives rise to the following 
expression for the effective gauge action,i2i2^ 

+ ^ J dtd'^ref'^^A^^d^A'^ 

+ dtd^e^'^Ald.Al (12) 



where 

Cc = 

Cs = 



247r2 
247r2 



|d3fcTr[Go^Go^Go^^""'- 



dk^ " 
-1 fl/^-i 



(13) 
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which are nothing but the charge (Cc) and spin (Cs) first 
Chern numbers of the colhnear magnetic ground state de- 
scribed by Go- The U(l) spin gauge field is defined as 
= Bfj^^z- It is straightforward to show that the fictitious 
spin gauge flux F^^, = dfj,A^ — d^A^ satisfies 



s d.Al^N ■[{d^N)x{d,N)]. (14) 



The effective gauge action in Eq. (fT2] i gives rise to the fol- 
lowing expressions for the charge current, 



dAC[r,t) ' 



271- 



47r 



A' 



(15) 



and the spin current, 



dAS{r,t) 
= ^''-'d.Af 



C, 



(16) 



In the absence of the spin gauge flux F^^^^, for example, 
for a homogeneous coUinear magnetic ground state, the above 
charge and spin currents describe the usual quantized charge 
and spin Hall effects driven by the external electromagnetic 
fields. The magnitude of the quantized charge (spin) current 
is determined by the charge (spin) Chern number. On the 
other hand, if the spin gauge flux can be generated by space- 
time dependent variation of the order parameter N through 
the relation in Eq. ( fT4] i. the adiabatic charge and spin Hall 
currents flow in the absence of the electromagnetic fields. In 
this case, interestingly, the magnitude of the quantized charge 
(spin) current is determined by the spin (charge) Chern num- 
ber in contrast to the case of the usual charge and spin Hall 
currents. 



B. Charge and spin pumping through a domain wall with 
spiral spins 

To demonstrate the idea of the spin gauge flux induced 
charge and spin responses, we consider the charge and spin 
polarizations induced by the adiabatic modulation of the spin 
direction in a domain wall of a CMI, which is described 
in Fig. EJa). An in-phase domain wall can be represented 
by iV(r,ry) in Eq. Q with 0(r) (j}{y) = f [-1 - 
2 tanh(y/^Dw)]- For convenience, we rearrange the compo- 
nents of N and consider N — Nx{^, t) given by 

Nx{r,t) = {cos 9 (t), sin 6 (t) cos (j){y), sin 0(t) sin(j){y)), 

(17) 

where the spin x-component is uniform. Traversing the do- 
main wall along the y direction, the spins rotate by 2tt within 
the yz-plane. External time dependent magnetic field H(t) — 
H{t)x coupled to the spins within the domain wall controls 



the magnitude of the uniform (x) components of the rotating 
spins. We assume that the length scale of the domain wall 
satisfies hv/Eg ^ ^dw ^ Ly. Here Ly represents the system 
size along the y-direction. Eg is the bulk energy gap, and v is 
the typical velocity scale contained in the bulk band structure 
near the chemical potential. Therefore the spin modulation 
near the domain wall is smooth enough but, overall, the spins 
are aligned collinearly throughout the entire system in which 
the spin z component Sz can be treated as a conserved quan- 
tity. 

Let us first consider the charge polarization induced by the 
spin gauge flux. From Eq. (fl4] l and Eq. (flST l. the adiabatic 
charge polarization current is given by 



on 

^ ^S,,:.N ■ [{dyN) X (dtN)]. 



(18) 



During the adiabatic evolution of 9{t) for i e [0, T], the spin 
electric field E^ = dyAf - dtAf^ =N ■ [{dyN) x (dtN)] 
produced by the spiral spins at the domain wall induces a 
charge current j^' , which leads to the average charge polar- 
ization (Px) = jj- J^" dy dtj^. In particular, if the unit 

vector N{y, t) fully covers the surface of the unit sphere in the 
(y, t) space during the adiabatic evolution, a quantized charge 
pumping is possible with the net charge transferred to the x 
direction given by 



Qpumped = j dt j dyj^ 



f de f ' d^N ■ [{d^N) X (dgN)], 



J7r/2 



5^skyrmion ; 



(19) 



where nskyrmion is the skyrmion number of the vector field 
N{y, t). Therefore the unit pumped charge is determined by 
the product of two topological invariants defined in momen- 
tum space (Cs) and position/time space (nskyrmion)- 

Moreover, since an electron carries both charge and spin 
quantum numbers, a charge current accompanies a spin cur- 
rent simultaneously. The spin current induced by the adi- 
abatic spin electric field can be understood in the follow- 
ing way. Since Sz is a good quantum number, a system 
with {C^,Ci) = {n^,Tii) can be considered as a super- 
position of two subsystems with (C^jCj.) ~ ("-t7 0) and 
(CV, Ci) — (0, ni), respectively. For given spin electric field 
Ey, the charge current induced in each subsystem should 

be — e^^^Ey and e-^^iJ^ , respectively, consider- 
ing the spin Chern number in each case. However, since the 
spin polarization directions are opposite in two subsystems, 
the total induced spin current is given by = ^^^^Ey + 

(-f =1 ^%+"^' ^^. Therefore the amount of spin 

transfer is determined by the charge Chern number Cc- The 
adiabatic spin current induced by the spin gauge flux can be 
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FIG. 5: (Color online) Charge and spin currents driven by the spin 
gauge flux which is produced by the slow variation of the uniform 
(x) component of the spins rotating in yz plane within a domain wall, 
(a) The rotating spins near the in-phase domain wall. Here the red 
rotating arrow describes the net ferromagnetic moment of a pair of 
the antiparallel spins within a unit cell. The external time dependent 
magnetic field H(f) — H{t)x couples to the net ferromagnetic mo- 
ment and controls the magnitude of the x-component of the rotating 
spins. In (b) and (c), the central square with graded colors in each 
panel, indicates the spiral domain wall producing the spin gauge flux 
on which the direction of the induced current is denoted with a red 
arrow. For a given spin gauge flux Ey, the direction of the adia- 
batic charge (spin) current is determined by the corresponding spin 
(charge) Chem number, (b) A TCMI with (Cs, Cc) = (1, -1). (c) 
A TCMI with (Cs, Cc) = (-1, -1). 



written as 

j;! = ^^'''^^ ■ [(^''^) ^ (^^^)]- ^20) 

The identical expression of the spin current can also be ob- 
tained from Eq. (fT4l i and Eq. (fTSI l. For an adiabatic cycle in 
which N{y,t) covers the full spherical solid angle, the net 
transferred spin along the a:;-direction is given by 



S'pumped = y dt J dyj^ 

<.-37r/2 

dd d<j>N ■[{d^N)x {dgN)], 

Jtt/2 

(21) 







Therefore the quantum of the pumped spin is determined by 
the product of the charge Chern number Cc and the skyrmion 
number riskyrmion- Interestingly, the charge and spin responses 
under the spin gauge fields are exactly dual to the correspond- 
ing responses under the electromagnetic fields. Given the spin 
gauge fields, the charge (spin) response is determined by the 
spin (charge) Chern number.— 



In Fig.|5](b) and (c), we show the adiabatic charge and spin 
currents induced by the adiabatic spin rotation at the spiral do- 
main wall for the two TCMIs with Cs = 1 and -1 in Fig.[T](b). 
Interestingly, while the charge cuiTents flow in the opposite di- 
rections because their Cs have opposite signs, the directions 
of the spin currents are the same in the two cases, consistent 
with the fact that Cc = — 1 in both cases. Adiabatic charge 
and spin currents for AHI can be obtained by adding the cor- 
responding currents for the two TCMIs, i.e. by the superposi- 
tion of Fig.|5](b) and (c). Since the charge currents flow in the 
opposite directions, the net charge polarization is zero. How- 
ever, the magnitude of the spin polarization is twice compared 
to the case of each TCMI. 



C. Charge and spin quantum numbers of a skyrmion 

It is worth noting that the last term of the effective gauge 
action in Eq. (fT2] i is nothing but the Hopf term, which deter- 
mines the spin and statistics of skyrmion solitons.^'' There- 
fore the topological invariants, which quantify the pumped 
spin and charge, also characterize the quantum numbers of 
skyrmion topological textures. The time components of the 
adiabatic currents in Eq. (fTSl l and Eq. (|20] | show the charge and 
spin densities associated with a skyrmion, which lead to the 
following expressions of the charge and spin quantum num- 
bers of a skyrmion, 

Qskyrmion = ^ / ^xdySi ■ {{p^N) X (dyN)], 



skyrmion 



Cc 



dxdyN ■ [{d^N) X {dyN)]. 



(22) 



Simple integration of the above equations shows that a 
skyrmion with the unit Pontryagin index on the TCMI is a 
spin- 1/2 fermion with the charge ±e. On the other hand, a 
skyrmion with the unit Pontryagin index on the AHI is a spin- 
1 boson with no electric charge. Therefore the skyrmion de- 
fects in TCMIs and AHIs carry nontrivial quantum numbers. 



IV. TOPOLOGICAL RESPONSES OF CHIRAL 
MAGNETIC INSULATORS DERIVED FROM SHI 

A. Effective action for the charge and spin responses 

The low energy Hamiltonian of a CMI which is derived 
from the SHI is given by 



i/eff= / d^rijHr,t)[{~td.,)T^v, + {-tdy) 



' l^sTz + t^TzVzOz + mKvTzNix, t) ■ a]4>{r, t). 



(23) 



In contrast to the case of CMI obtained from the AHI, the lo- 
cal SU(2) spin rotation of the staggered spin order parameter 
N{r,t) cannot make the Hamiltonian to be diagonal in spin 



space because of the simultaneous rotation of the spin depen- 
dent <2 term representing the strength of the spin-orbit cou- 
pHng. To obtain the effective gauge action through the gra- 
dient expansion, it is necessary to transfer the position/time 
dependance of the Hamihonian to the gauge fields via local 
unitary transformations, which can be achieved only in the 
limit of jmApl ^ 1^21 or IttiafI ^ 1^21 corresponding to SHI 
or SDW phases in Fig.[3](c), respectively. Therefore we focus 
on the adiabatic responses of SHI and SDW phases in the fol- 
lowing discussion. Moreover, since the spin is not conserved 
in general due to the spin-orbit coupling, we only consider the 
charge polarization currents. 

We first divide the Hamiltonian into two parts in such a way 
as i/eff = ^^eff,+ + ^^eff,- in which 7Jeff,+ (Heff,.) corresponds 
to = 1 (i^z = — 1) describing the low energy particles near 
the node at k = Ki (k = — Ki). Explicitly, H^ff ± are given 
by 

H,s,± = J d\i:i{r,t)[±{-td.,)T, + {-idy)Ty 

+ fJ-sTz +T^{±t2a:, +mAFN{r,t) ■ a}]ilj±{r,t). 

(24) 

To make the Hamiltonian to be diagonal in spin space at each 
node separately, we define the collective spin degrees of free- 
dom N± which combine N ■ a and t2crz in the following way, 

± ^20-2 + rnApN{r, t) ■ a, 

= mp^N^a^ + friAFNyay + {mp^N^ ± t2)az, 

+ tl± 2mAFt2N,N±{r, t) ■ a, (25) 

where the new unit vector N± satisfies 

7V± • [(a^iV±) X (9,7V±)] 



(a) 



[m2p + i2±2mAFi2A^2]3/2 



(26) 



Now we consider the local unitary rotations \J± of the 
unit vectors N± satisfying C/j_((T • N±)\J± = a^- Defining 
■\\j± = C/iV^, the effective action can be summarized as that 
in Eq. ( fTOb but with G^^ replaced by 

Gg ^(fc) =W - (±)fc^Ta. - kyTy - ^sTz 



+ tl±2mAFt2NzTzaz, (27) 



and ia • if^ = iUl{df,U+) or iUl{d^U^). In Eq. (HTJ, Gq 
depends on the position/time coordinates due to the presence 
of Nz term. Therefore when the variation of Gq in the posi- 
tion/time spaces becomes important, the gradient of Gq can 
also make a nontrivial contribution to the charge polarization 
as the gradient of the charge/spin gauge fields does. Therefore 
the effective gauge action approach in the previous section is 
valid only when jmApl ^ 1^21 or jmApl ^ 1^21 is satisfied, 
in which has a negligible contribution to Gq. However, 
regarding the charge polarization currents shown in Eq. ( fTST i 
and Eq. ( fTSl l. Gq appears only in the definition of the Chern 
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FIG. 6: (Color online) Spin structures at anti-phase domain walls 
for mAF = tuafZ- Here the arrows indicate the staggered spins of 
the antiferromagnet (or ferrimagnet). Spins rotate with respect to an 
axis, X or y, which is normal to the collinear spin direction z. (a) The 
Neel domain wall in which spins rotate around x axis, normal to the 
spin modulation direction y. The spin structure can be described by 
using iV^* in Eq. i30i . (b) The Bloch domain wall in which spins 
rotate around y axis, parallel to the spin modulation direction. The 
spin structure can be described by using TV^^ in Eq. 



number, which is topologically invariant. Therefore as long as 
Nz does not induce the change of the Chern numbers via band 
gap-closings, dependence of Go does not affect the charge 
polarization current obtained by considering the gradients of 
gauge fields. 

Through the gradient expansion of gauge fields following 
the same procedure as in Sec. lIII Al the adiabatic charge cur- 
rent can be obtained as 



Cs3 
Stt 



■e^,xN± ■ [{d,N±) X {dxN±)], (28) 



where Cs.± — Tl for both jmApl > 11^21 + ImsII ™d 
I'tiafI < 11^2! — ImsII cases corresponding to SDW and SHI 
phases, respectively. The total charge polarization of the sys- 
tem can be obtained by adding the adiabatic currents from the 
two nodes. 



- jC 



(29) 



The adiabatic charge polarization (and also the skyrmion 
charge) for arbitrary ItoafI, especially for the TCMI phase 
correspondingto ||f2| — ImsII < \^af\ < ||i2| + lMslMs inves- 
tigated in Sec.[V]by considering the systematic semiclassical 
expansion for adiabatic topological responses of inhomoge- 
neous crystals. 



B. Charge polarization at a domain wall with spiral spins 

We consider the charge polarization or pumping through a 
domain wall with spiral spins. Because of the spin anisotropy 
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due to the spin-orbit interaction t2 = t2Z, the relative orien- 
tations of mAF, t2, and the rotation axis of domain wall spins 
are crucial ingredients determining the charge polarization. To 
describe the domain wall spins, we consider the following rep- 
resentations of inhomogeneous spins, 

7V^*(r,t) = (cos6l(i),sin6'(t) cos (/)(y), sin 6l(t) sin (/)(y)), 

N^^{Y,t) = {sme{t) sin(j){y),cose(t),aine{t)cos(j){y)), 

N^^{r,t) = {sme{t) cos (j){y), sme{t) sin (l){y), cos e{t)). 

(30) 

In iV°*, the spin a-component is cos 9{t) which is uniform 
independent of y. When we describe a domain wall with 
iV°*, a indicates the rotation axis of the domain wall spins. 
The deviation of 9 from ^ describes the adiabatic change of 
the relative magnitude of the uniform and rotating components 
of the domain wall spins. 

In Fig. |6] anti-phase domain walls with rotating spins are 
described for mAp = rriAFZ- Here with the boundary 
conditions of (f){y — oo) = ^ and (j){y +oo) = de- 
scribes the Neel domain wall in which the spins rotate around 
X-axis perpendicular to the spin modulation direction y as 
shown in Fig.|6](a). On the other hand, TVy* with the bound- 
ary conditions of (f>{y —oo) — and +oo) = tt, 
describes the Bloch domain wall in which the spins rotate 
around the spin modulation direction y as shown in Fig.|6](b). 
For these two configurations we compute the partial polariza- 
tion that is defined as 

(31) 

Using Eq. ( |26] l, Eq. (|28] l, and Eq. ( |29] l, it can be easily shown 
that ^ for both 7V^* and iV^?* are zero for any 9. There- 
fore no charge polarization occurs at the domain walls with 
rotating spins when mAp || t2. 

Similarly, the partial polarization induced by the rotat- 
ing spins at the domain walls can be considered for mAp -L t2 . 
If mAP II X, iVf* (iV^*) describes the Neel (Bloch) domain 
wall with appropriate boundary conditions. It can be shown 
that ^ ^ in general for the Neel domain wall (iV^*). 
On the other hand, in the case of the Bloch domain wall de- 
scribed by TVy*, has opposite signs for two inequiva- 
lent domain wall configurations in which (f>{y) E [— ^, ^] or 
(f>{y) G respectively. Therefore the average polariza- 

tion for this case is expected to be zero when two anti-phase 
domain wall configurations are equally distributed in material. 
If mAP II y, both iV^* and iVf* describe the Neel domain 
walls. For the Neel domain wall described by iVf ^ 

in general. However, in the case of the domain wall described 
by N^, has opposite signs for two inequivalent domain 
wall configurations in which e [0, vr] or (j){y) £ [vr, 2tt], 
respectively. The average polarization for this case is also ex- 
pected to be zero. 

Therefore the partial charge polarization can be 

nonzero only for mAp -L t2 at the Neel domain walls de- 
scribed by TV^*. from the Neel domain wall is shown 
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FIG. 7: (Color online) Charge polarization at the Neel domain wall 
described by iVf* = (sin 6l(t) cos <?!)(?/), sin 6l(t) sin (^(y), cos 6l(t)) 
for mAF II x. Here the boundary condition is given by </>(y — >■ 
— oo) = and (j){y — >■ oo) = n. (a) Spin configuration near the 
Neel domain wall, (b) as a function of 9. 

in Fig. |2] Interestingly, if we assume that the spin rotation 
occurs over the whole lattice system, 7V^* (iVy*) describes 
the transverse (longitudinal) conical magnetic ground state. — 
Therefore the charge polarization occurs only for the trans- 
verse conical magnetic phases in the direction perpendicular 
to the spatial modulation, consistent with the general theory 
for the charge polarization in spiral magnets. ^''•^^ 

Now let us discuss about the quantized charge pumping 
(also the charge quantum number of a skyrmion) consider- 
ing the variations of e [0, 27r] and 9{t) E [0,7r]. At first, 
for iV°*, the local adiabatic charge currents j^^{r,t) — 
j±. ^{9, (p) can be obtained, which satisfy the following rela- 
tions, 

£j9,4,) = -j^Jn-9,ct>). (32) 

Because of this relation, the total charge current j^{9, (j)) — 
^{6, (j)) + ^{9, (j)) also satisfy the following condition, 

jC{9,ct,) = -j<^{^-9,4,). (33) 

Therefore if the uniform component Nz — cos 9{t) of iV(r, t) 
changes over the full polar angle 9{t) G [0, tt], the net trans- 
ferred charge due to the adiabatic charge current j^' is zero, 
i.e., 

Qpumped ^ dtj dyj^ = 0. (34) 

On the other hand, for N^{r,t) and iV^*(r,t), it can be 
shown that the adiabatic charge polarization is always zero 
for any given 9{t) after the integration over 0, because of the 
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following condition satisfied by 



Jx,- 



,(/) + 7r). 



(35) 



Therefore the quantized charge pumping is impossible for any 
orientation of biaf- 

The fact that the net transferred charge is zero can be un- 
derstood in terms of the spin gauge flux induced spin Hall 
effect. If I TO AF I ^ 1^2! where /V+ and iV_ are almost paral- 
lel, the charge current density is proportional to C5.+ + Cs,- 
as shown in Eq. ( |28] | and Eq. ( |29] |. Since Cs.+ = —1 
and Cs,- = 1 in this limit, the total charge current is zero 
for given spin gauge flux. The opposite sign of Cs,± is 
the direct consequence of the fact that the Dirac particles at 
the two nodal points have the opposite winding directions, 
which is reflected in the opposite sign of the dispersions 
in Eq. dm.^"'^' However, when |toaf| is finite, and iV_ 
are tilted with a small relative angle due to the spin-orbit cou- 
pling. In this case, since the spin gauge fluxes generated by 
N± are different, the local polarization current can be nonzero 
although Cs,+ — ^1 ™d Cs,- — 1 ^re still satisfied. There- 
fore the finite partial polarization is essentially the result 
of the spin anisotropy due to the spin-orbit coupling, not of a 
topological origin. On the other hand, if ItoafI ^ 1^21 cor- 
responding to SHI phase, the solid angle subtended by N± 
is very tiny even though N covers the full spherical solid an- 
gle, which leads to vanishingly small charge polarization. The 
quantized pumping is also impossible in SHI case. 



C. Skyrmion charge 

The skyrmion defects of the colUnear magnetic phases sat- 
isfying mAF II t2 or mAF -L h can also be described by using 
the spin configurations z * l30l l. For the descrip- 

tion of a skyrmion, we replace 9{t) by 0{\r\) satisfying the 
boundary conditions 6'(|r| = 0) = tt and ^(Irl — > 00) = 
and adopt (/){r) = tan^^(^) with r ~ {x,y). In the case of 

mAF II t2, a skyrmion defect can be described by using iV^*. 
As Eq. dSTT l. Eq. ( l33T l and Fig. |7](b) show, the adiabatically 
induced charge of a skyrmion, which is equivalent to the par- 
tial polarization of a domain wall, is anti-symmetrically dis- 
tributed with respect to {6 ~ tt/2). Therefore the integration 
of the induced charge adds up to zero, which shows that the 
skyrmion is not charged. On the other hand, when mAF -L t2, 
a skyrmion can be described by or Ny^ . In this case, as 
shown in Eq. ( |35] l, the induced charge is zero for any 9 value. 
Therefore the skyrmion is not charged. 

To sum up, the skyrmions in CMls derived from SHI do 
not carry the charge quantum number at least for |toaf| < 
11^21 — |a*s|| or ItoafI > 11^21 + \fJ-s\\, which corresponds to 
SHI or SDW phases, respectively. 



D. Meron charge 

The anti-symmetric distribution of the induced charge with 
respect to {6 — tt/2), shown in Fig.|2](b), implies that a meron- 



type defect in which 9 changes over 9 € [0, 7r/2] can be 
charged although a skyrmion is chargeless. For mAF || ^2, a 
meron can be described by iV°* with the boundary condition 
given by 6'(|r| = 0) = tt/2 and 6'(|r| 00) = 0. Explic- 
itly, from the adiabatic charge current in Eq. ( |28] |. the meron 
chai-ge Qmeron is givcn by 



1 



t2 



'AF 



V 



TOaf + 1 



2tl ' 



(36) 



for SHI phase satisfying t2 ^ toaf > 0. Similarly for SDW 
phase which satisfy toaf ^ ^2 > 0, 

t2 t2 



Qn 



Toaf 



(37) 



Therefore a meron defect has a finite charge, however, which 
is not quantized. 

A meron defect is especially important for the spin order- 
ing with easy-plane anisotropy. Considering a planar antifer- 
romagnetic spin ordering mAF = {'m-x,my, 0), the low energy 
Hamiltonian can be written as 

i?eff,± = j (fr'ijj!^{r,t)[±{-idx)Tx + [-idy)Ty + ns^z 

+ Tz{±t2(7z + rrixGr^ + my(7y)]^±[Y,t), (38) 

where ± indicates the two nodes at k = ±Ki. After a trans- 
formation of V'- Ty(Tzip-, the full Hamiltonian iJetf = 
Heff,+ + i?eff,- can be written as 

Heff = J cPr^Hr, t)[{--idx)Tx + {~idy)Ty + ^isTyV^ 

+ Tz{t2<Jz + mx(7x + myay}]ip{r,t), (39) 

where 1]:^ ~ (?/'^,^/'^). It is interesting to notice that this 
Hamiltonian has similar structure as the effective Hamiltonian 
of CMI derived from AHI in Eq. Q. Now we define a three 
component unit vector TV — / 2 ^ 2 2 ("^a;, n^y, ^2) and 

compute the charge current induced by the adiabatic variation 
of N using the formulation in Sec. IIII A I A meron can be rep- 
resented by 

^mei-on — (cos a(r) COS 0, COS a(r) sin 0, sin Q!(r) ) , (40) 

which is under the constraint of a(r = 0) = ^ and sin[Q;(r 
oo)] = *^ Here to^ = to^ + ml, ^(r) = tan~-^(-), 

and |r| = r = + y^. If |/is| > \/ m? + t\, the induced 
charge is zero. On the other hand, when |/is| < \/ m? + t\ 

ron del 

[{dxN. 



the total charge induced by a meron defect is given by 



4tt 



dxdyN„ 



t2 



- 1, 



X ( Oy Af^meron ) ] ; 

(41) 



which shows that a meron defect is charged but the net charge 

2 

is not quantized. Qmei-on ~ ^p- for ^2 > while Qmemn ap- 
proaches -1 for m ^ t2- The meron charge of the easy-plane 
Neel phase derived from SHI is also discussed by the recent 
work in Ref. 32, in which the same conclusion is obtained. 
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V. SECOND CHERN NUMBER AND TOPOLOGICAL 
CONTRIBUTION TO THE INHOMOGENEITY INDUCED 
ADIABATIC CURRENTS 

A. Semiclassical gradient expansion approacli for adiabatic 
currents 

In this section, we demonstrate that the charge and spin po- 
larizations induced by the adiabatic variation of the inhomo- 
geneous spin order parameter A^(r, t) can be understood as a 
generic topological response of inhomogeneous systems un- 
der adiabatic temporal variations. In particular, it is shown 
that the topological nature of the adiabatic current of inhomo- 
geneous systems is endowed with the second Chern number. It 
turns out that the existence of the second Chern number is the 
topological origin of the fact that the spin and charge polariza- 
tion quantums induced by the inhomogeneous spin order pa- 
rameter N{r, t) are determined by the product of two topolog- 
ical invariants, which are defined in the position and momen- 
tum spaces, respectively, as shown in Eq. (T% and Eq. ( 1211 1. 

For simplicity, we assume that the system has translational 
invariance along the x-direction and the spatial inhomogene- 
ity is developed along the y-direction. The total charge current 
along the x-direction, is given by 



Jx = -— — iXrlnG 



-le 



27r 



dy I dttr{t,y\G{K)dk^G-\K)\t,y), 

(42) 



where the symbol tr indicates the trace over the discrete in- 
dices while Tr includes the trace over discrete indices and in- 
tegration over the momentum and the position y and time 
t. Since the electrons are coupled to the time-dependent in- 
homogeneous order parameter N{y, t), the Green's function, 
which is diagonal in k^, has off-diagonal components in the 
{y,t) space, i.e., (yi, ti|G(fc:r)|y2, ^2) = G{k^;yi,ti;y2,t2). 

If N{y, t) changes slowly in the position and time spaces, 
namely, if Ly ^ ^/Eg and 1/T ^ Eg are satisfied, the 
position and time dependence of the Green's function G can 
be treated by the gradient expansion method. Here Ly (T) 
is the length (time) scale over which N{y,t) varies and Eg 
is the bulk energy gap. To perform the systematic gradi- 
ent expansion, we first consider the Wigner transformation of 
G{kx;yi, ti;y2, ^2), i-e- the Fourier transformation with re- 
spect to the relative coordinates yi — j/2 and ti — <2, which 
can be written as"'-'"-'^ 

p^., , yi + y2 ^l+<2^ 

G[kx]ky, ;w, — - — ) 



X G{k^\y2,ti;y2,t2). 



(43) 



As shown in detail in the Appendix lAl the key ingredient of 
the gradient expansion is that the Wigner transformed Green's 
function G can be expanded order by order in powers of 



the spatial gradient of the semi-classical Green's function 
Q{kx;ky,y\uj,t) ^ [uj - 'H^^{k^;ky,y]U},t)Y'^?^'^^ Iniht 
semi-classical Hamiltonian Tiscikx', ky, y; oj, t), the conjugate 
variables {y, ky) and are treated as independent real 

numbers {y, ky) and (i, w). Here ky — and Cj — i^. 

Therefore the gradient expansion is equivalent to the semi- 
classical expansion. 

The leading order term of the adiabatic current induced by 
the spatial inhomogeneity can be written as 



Jx — <-^2:, topological ^ <-^a:,peitmb; 



(44) 



in which 

J X , topological 



du 



dydt 



60 J 27r (27r)2 



(45) 



where the indices {I'l, 1^2, 1^3, 1^4, 1^5) run over {u, k^, ky,y,t). 
Here indicates the 2-torus made of two dimensional Bril- 
louin zone. If the {y, t) dependence of the Hamiltonian oc- 
curs only through the coupling to iV(y, t), J^. topological can be 
rewritten using the spherical coordinates (6*, 0) for the unit 
vector N{y, t) = N{6{y, t), (t){y, t)) in the following way. 



J. 



X, topological 



le 
~ 60 

X £ 



did 

2^ 



T2 (27r)- js 



=eC, 



(2) 
Chern 



(46) 



run 



over 



where the indices (/ii, ^2, M3; M4j A^s) 
{uj,kx,ky,0,(f)). When N{y,t) fully covers the unit 
sphere during the adiabatic variation, the total adiabatic 

current is quantized with the magnitude of the quantum 

(2) 

determined by the second Chern number C^jjcm- other 
non-topological terms, whose contribution would vanish for 
cyclic evolutions of N, are included in Jj^^pertmb- The explicit 
expressions of Ja^^penurb are given in Eq. (IA28b . Eq. ( IA29b . 
and Eq. (|A30] l. 

The charge density induced by the spatial inhomogeneity 
can also be derived in a similar way from the semiclassical 
gradient expansion approach. Starting from the expression of 
the integrated charge density which is given by 



Jo = -— -iTrlnG 

= -ie — I dx I dyti{x,y\G{uj)du;G~^{uj)\x,y), 



(47) 



the inhomogeneity induced charge can be obtained from the 
expansion of the Wigner transformed Green's function in the 
powers of the spatial gradients and dy. The lowest order 
contribution to the adiabatic charge induced by a skyrmion 
defect appears in the second order terms which contains dx 
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and dy at the same time. Explicitly, it is given by 



where = + /| 



Then 



Qskyrmion 



le 
60 



2^ 



X £ 



mfJ.2fJ.3fJ.ifJ-'. 



T2 (27r)2 jg: 

,,tr[g(a,,,g-i: 



(48) 



(51) 



Using ^ from Eq. ( |49] >, it can be shown that 

9 



■EabcU^acU. (52) 



Taking into account the fact that is spin-diagonal, it can be 
shown that 



eabcNcin[a:,gD (Sfc^ gD^)gD {dti.^gjj^)gD {du^gi:,^)] : 

(53) 



dy dt 



where the indices (/ii, /i2, /is, /i4, /is) run over 
{uj,kx,ky,d,(j)). Here 9 and (p represents the spher- 
ical coordinates for the three component unit vector 
N{x,y) ~ N{9{x,y),(j){x,y)) describing a skyrmion con- 
figuration. The same expression for the skyrmion induced 
charge was obtained in the recent work by Santos et al. 
in Ref. [33 It is interesting to notice that Qskyimion can be 
obtained simply by replacing the time coordinate t by the 
spatial coordinate x from the expression of the total adiabatic 
charge current J^, topological in Eq. (gill. 

Now we show that the appearance of the second Chern 
number in J^. topological is the topological origin of the fact that 
the quantum of the pumped charge in CMls is given by the 
product of two topological invariants as shown in Eq. ( fT9] l. 
Let us consider the following semi-classical Green's function 
g which describes the low energy excitations of CMIs^^ 

g-^=u- /i(k)r, - /2(k)Tj, - mT,N{y, t) ■ a, (49) 

where r^.y.z and ax,y,z are Pauli matrices and /i,2(k) are 
functions which depend only on the momenta k. Let us con- 
sider the property of the following term, 

= tr[g{dN^g-^)g{dN,g-^)g{dk,g-^)g{dk2g~^)gid^g~% 

(50) 

where Na is a component of the unit vector N = 
{N.j:, Ny, Nz). We introduce a local unitary transformation 
U that satisfies ■ N)U — cfz- The semi-classical 

Green's function in the rotated frame, which is defined as 
go = WgU, is diagonal in the spin a space. After insert- 
ing the identity U'^U = 1 between g and g^^, lNaNtkik2Uj 
transforms as 

iN^N,k,k2u = ti[u^g{dN^g-^)g{dN,g-^)u 

X gD{dk,gD')gD{dk2gD')gD{dugD')] 



X tr[<JzgD {dk, gD')gD {dk2 gD^)gD (a^e^' )] . 



(54) 



Interestingly, the inhomogeneous part that depends on (y, t) 
and the homogeneous part that is written in term of go in 
the {uj,kx,ky) space are completely separated. Applying this 
result to Eq. ( |45] l. J^r, topological can be written as 

J.,topoiogicai = ^ / dydtN ■ [{dyN) X (dtN)] , (55) 

where the spin Chern number C5 is defined in terms of gjj 
using Eq. ( flj] ) after replacing Go by tJc- It is to be noted 
that the total charge current Ja;. topological induced by the adia- 
batic change of the spatial inhomogeneity is equal to the net 
transferred charge, which is shown in Eq. iT% . The existence 
of the topological contribution to the inhomogeneity induced 
charge current characterized by the second Chern number is 
the origin of the quantized charge pumping with the quantum 
of the pumped charge expressed as the product of two distinct 
topological invariants. 

Following the similar line of reasoning, the inhomogeneity 
induced total spin current can be obtained in the following 
way. 



J: 



X, topological 



X 



le 
'60 



duj r d^k 

2^ Jt2 (27r)2 



dydt 



,.,tr[(-a.iV(y,0)e(9.,a- 



■g{d.,g'% 

(56) 



For the semi-classical Green's function in Eq. ( |49l ) it is 
straightforward to show that. 



J: 



a;, topological 



/ dydtN ■ [{dyN) X (dtN)] , (57) 



where Cc is the charge Chern number defined in terms of g. 
It is interesting to notice the equivalence of topological ^bove 
and the total transferred spin Snansfer in Eq- <E}- 

Recently, the topological contribution to the inhomogene- 
ity induced adiabatic charge polarization was investigated by 
considering the Berry phase approach to the charge polariza- 
tioni^^ Using the semiclassical wave packet dynamics for- 
malism, Xiao et al.," '^ have derived the general expression 
of the adiabatic polarization induced by spatial inhomogene- 
ity. For a 2D multi-band system in which the smooth spatial 
modulation exists along the (3 direction (/? = x,y) with the 
modulation period of L^, the net charge polarization along 
the a direction perpendicular to the modulation direction (3 is 
given byiiril. 
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de 




FIG. 8: (Color online) Inhomogeneity induced adiabatic polariza- 
tion (Px) for a CMI derived from AHI. A spiral spin ordering with 
the modulation period of Ly along the y-direction is assumed. 
is plotted in the unit of e using Eq. l l58l l. For t2=2, /is=l, is 
computed for mAF=-2, 0.5, 2 and 4 which correspond to a TCMI 
(Cs = -1), an AHI (Cs = 0), another TCMI {Cs = 1) and a 
SDW (Cs = 0), respectively. Notice that the polarization is solely 
determined by the spin Chem number Cs- 




FIG. 9: (Color online) Inhomogeneity induced adiabatic polarization 
(Px) for a CMI derived from SHI. A spiral spin ordering with the 
modulation period of Ly along the j/-direction is assumed. is 
plotted in the unit of e using Eq. l l58t . For t2=2, /is=l, is com- 
puted for mAF=0.5, 2 and 4, which are relevant to SHI, TCMI, SDW 
phases in Fig.[3](c), respectively, in the collinear limit. 6 = dc and 
9 — {tv — 6c) indicate the locations where the singularity in 
occurs for mA¥=2. 



where (r, t) is the inhomogeneity induced adiabatic current 
as the parameter \{t) evolves from A(0) to A(r). Here the in- 
dices z, J, fc, I run over [ka ,ki3,ri3,X) and the Berry curvature 
Fij satisfies Fij = diAj — djAi — i[Ai, Aj] with the Berry 
connection A^"^ = {^f^\idj\^,y). Here \^^) is the periodic 
part of the Bloch wave function for the occupied band i^. The 

equivalence of the Pi^' in Eq. ( |58] l and Jrc, topological in Eq. ( |45] l 
can be proved by considering the topological invariance of the 
second Chern number when the order parameter N{y,t) spans 
a closed manifold. The details of the proof is given in Ref. 
Therefore the inhomogeneity induced charge polarization can 
be computed, for a general (lattice) Hamiltonian, using either 
Eq. (US) or Eq. 



B. Adiabatic polarization and skyrmion charge for CMI 
derived from AHI 

To demonstrate the adiabatic polarization induced by inho- 
mogeneous ferrimagnetic ground states, we consider a spiral 
(or conical) spin ordering pattern with the modulation period 
Ly = 2tt/Q. The periodic spin modulation along the y di- 
rection can also be represented by use of z '^l- •EHJ 
taking = Qy. The spatial modulation of the spin order- 
ing is described in Fig.|4] The polarization of the full lattice 
Hamiltonian iJfuu in Eq. Q is numerically computed from 



the second Chern form given in Eq. ( fSST l. Since the modu- 
lated spin ordering spreads over the whole system, Sz is not 
conserved any more. Therefore we focus on the charge polar- 
ization induced by the inhomogeneous spin ordering. 



Let us first consider the CMI derived from AHI. In Fig. [8] 
we plot the derivative of average polarization per unit 
modulation period Ly as a function of the polar angle 9{\) 
which plays the role of the adiabatic parameter. For t2=2, 
/is=l, the polarization is computed for mAF=-2, 0.5, 2 and 4 
which, in the collinear limit, correspond to a TCMI (Cs — 
-1), an AHI (Cs = 0), another TCMI (Cs = 1) and a SDW 
(Cs = 0), respectively, following the phase diagram in Fig.[T] 
(b). It is worth to be noted that the magnitude and sign of the 
polarization are solely determined by the spin Chern number 
Cs- Moreover, as 9 evolves from to tt, the average polar- 
ization per unit modulation period Ly is quantized satisfying 
Jq dO^^ = Cs- Therefore quantized charge pumping is pos- 
sible in CMI derived from AHI, consistent with the conclusion 
in Sec. |III] At the same time, this means that a skyrmion de- 
fect made of the staggered spin N has the charge quantum 
number Cs- The polarization (or skyrmion charge) is also 
computed including the ferromagnetic component of the spin 
order parameter rrip considering the general ferrimagnetic or- 
dering with a net ferromagnetic moment. We have confirmed 
that mp does not affect the polarization as long as the bulk 
band gap is not closed by rnp- 
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C. Adiabatic polarization and skyrmion charge for CMI 
derived from SHI 



In the case of the CMI derived from SHI, depends on 
the orientation of the uniform component cos 9 of the rotat- 
ing spins relative to the anisotropy direction t2 inherent to the 
SHI. In the case of the spiral ordering described by or 
TVy* in Eq. ( l30b . in which the uniform component aligns per- 
pendicular to the spin anisotropy direction, — always. 
Therefore the net pumped charge is zero in this case. 

On the other hand, for the spiral ordering described by 
in Eq. ( l30t , which corresponds to the transverse conical phase 
with the uniform spin component parallel to t2, is plotted 
in Fig.|9] For i2=2, /i5=l, the polarization is computed for 
?7iAF=0.5, 2 and 4, which are relevant to SHI, TCMI, SDW 
phases in Fig.[3](c), respectively, in the coUinear limit, ^jg- is 
odd function with respect to 6* = ^, which is consistent with 
Eq. ( |33] |. Because of this, as evolves from to tt, the net 
polarization is zero, i.e., /J^ dO^^ ~ 0. Therefore the quan- 
tized charge pumping is impossible via the adiabatic rotation 
of the uniform component of the spiral spins in this case. At 
the same time, it means that the skyrmion defect has no elec- 
tric charge independent of |toaf| in the case of CMI derived 
from SHI. 

The CMIs with m\p=2 satisfy the condition of |i2 — /-isl < 
I'tiafI < 1^2 + /^s|> in which the adiabatic polarization could 
not be properly described by the effective gauge field ap- 
proach used in Sec. |IV] In this case, shows strong en- 
hancement near 9 — 6c and n ~ 9c as shown in Fig. |9] At 
the singular points of 9 = 9c and tt — 9c, is not well- 
defined. Comparing the magnitude of near 9 = 9c or 
9 = TT — 9c, the partial polarization for m\p=2 is an order of 
magnitude larger than those for toaf = 0.5 or toaf = 4. The 
origin of this enhancement can be understood in the follow- 
ing way. The semiclassical Hamiltonian describing the low 
energy electrons near k = ±Ki in CMI derived from SHI is 
given by 

(59) 

Using the semiclassical Green's function Q^^ = uj — Use, the 
local charge Chern number can be defined in the following 
way. 



Cc{r,t) 



247r2 J dk^ dK dkx " 

=C+(r,t) + C_(r,0, (60) 



in which 



C±{r,t) = (t)^ sign^^s + \/ + *i ± 2mAFi2iV, 



sign! A^s - V "^Af + ^2 ± 2mAFt2N2 



(a) 




(b) 



(c) 



Cs=0 






FIG. 10: (Color online) The distribution of local charge Chern num- 
bers within a skyrmion in the case of (t2 — f-is) < rriAF < {t2 + /is) 
and t2 > i-is > 0. Here the skyrmion configurations for different 
orientations of hiaf relative to t2 — t2Z are considered separately. 
The largest dotted circle describes the size of a skyrmion. Red solid 
lines indicate the phase boundary where the charge Chern number 
changes due to gap-closing, (a) For hiaf = rriAFZ (t2 || hiaf). (b) 
ForniAF = niApx (t2 -L hiaf). (c) ForniAF = m^vy (t2 1- hiaf). 



In contrast to the cases of |toaf| > 1^2 + Msl or IttiafI < 1^2 — 
/i^l in which Cc(r, t) — Q for any r and t, when \t2 — /is| < 
I'tiafI < 1^2 + ^^s\, Cc'{r,t) changes discontinuously in the 
(r, t) space because the energy spectrum of the semiclassical 
Hamiltonian Hsc develops band touching points. At 6* = 6'c 
and n ~ 9c where develops singularities. Use possesses 
gapless points at k = Ki or k = — Ki, which is the origin of 
the strong enhancement of However, since the systems 
becomes gapless at this point, leakage currents will flow in 
reality, which disturbs the charge polarization. 

The properties of skyrmions in CMI satisfying \t2 ~ fis\ < 
|"1af| < 1^2 + /^s I can also be extracted from the spatial 
distribution of the local charge Chern number Cc(r, t) = 
Cc{x,y). We consider the following simple skyrmion con- 
figuration. 



Nix,y) = i- 



2x 



2?/ 



y 



1 



+ 1' 



1' 



(62) 



(61) 



which corresponds to the case of hiaf || t2 for t2 = t2Z. 
For (^2 — Ms) < ^^AF < (^2 + ^J-s), the local charge Chern 
number has spatial variation whose distribution is shown in 
Fig. Uni (a)- The red solid circle connects the points where 
the bulk gap vanishes locally and the Chern number jumps. 
Due to the occurrence of the gap-closing, the skyrmion defect 
cannot be constructed from the adiabatic deformation of the 
spin ordering directions and it is not energetically favorable to 
make skyrmions in this phase. 
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Because of the spin anisotropy induced by t2, the distribu- 
tion of the local charge Chern numbers also depends on the 
orientation of mAp- The local Chern number distribution for 
t2 -L mAF, which can be described by permuting the compo- 
nents of N{x, y) in Eq. ( |62] ). are shown in Fig.[TO](b) and (c). 
The rotation of mAp relative to t2 changes the location of the 
circles. The size of the solid circle depends on the magnitude 
of TTiAF- The topological equivalence of the distribution of the 
local Chern numbers in Fig. [10] (a), (b), and (c) can be seen 
after compactifying the boundary of the skyrmion (the dotted 
circle) into a single point. Therefore a skyrmion configuration 
always accompanies the band gap-closing. 



VI. CONCLUSION 

In this work, we have investigated the adiabatic charge and 
spin polarizations in chiral ferrimagnetic insulators through 
the smooth tilting of the inhomogeneous spin structures. 
There are two essential ingredients for the adiabatic polariza- 
tions. One is the bulk topological property and the other is 
the spin anisotropy. In the case of CMI derived from AHl, 
the charge and spin polarizations are solely determined by the 
bulk topological properties represented by the charge and spin 
Chern numbers. The charge and spin Hall effects induced by 
the spin gauge flux, which can be interpreted as dual responses 
compared to the usual Hall effects, drive the charge and spin 
polarizations. Because of its topological nature, quantized 
pumping is possible and a skyrmion carries quantum numbers 
in this CMI. In particular, when the system is lack of the spa- 
tial inversion symmetry, a new CMI phase dubbed the topo- 
logical chiral magnetic insulator (TCMl) exists. In TCMls, 
the quantized charge and spin pumpings are possible and the 
skyrmion texture is a spin- 1/2 fermion with the charge e. 

On the other hand, spin anisotropy plays the crucial role in 
the case of the charge polarization in CMls derived from SHI. 
Because of the spin anisotropy induced by the spin-orbit cou- 
pling, the charge polarization occurs even in the topologically 
trivial SDW phases. The charge polarization is also strongly 
affected by the bulk topological property. When the magni- 
tude of the spin ordering becomes comparable to the spin-orbit 
coupling strength which corresponds to TCMI phase in the 
collinear limit, the local charge Chern number changes dis- 
continuously in the position/time spaces. Approaching the lo- 
cal phase boundary where the discontinuity in the local charge 
Chern number occurs, the charge polarization shows a strong 
enhancement due to the existence of gap-closing points. In 
addition, in the case of a meron type defect, it can be charged 
in contrast to the case of skyrmions. In particular, the meron 
of the easy-plane Neel state carries finite charge but the net 
charge is not quantized. 

The quantum number of a skyrmion, which exists in CMI 
derived from TBI (both AHI and SHI) can be summarized 
in the following way. If the collinear magnetic ground state 
has finite Chern numbers independent of the spatial orienta- 
tion of the spins and supports a stable skyrmion defect, the 
skyrmion carries a nonzero quantum number In the case of 
CMI derived from AHI, the relation of the skyrmion quantum 



number with Chern numbers of the collinear magnetic ground 
state is obvious. The charge quantum number of a skyrmion is 
given by the spin Chern number. On the other hand, the spin 
quantum number and quantum statistics of a skyrmion are de- 
termined by the charge Chern number. In the case of CMI 
derived from SHI, since the spin Chern number is not well- 
defined, the topological property of CMI is determined by the 
charge Chern number. In both SDW and SHI phases, since the 
charge Chern number is zero, there is no topological contribu- 
tion to the skyrmion quantum numbers. On the other hand, in 
the case of TCMI, although it has a finite charge Chern num- 
ber, a skyrmion is not a stable object because it always accom- 
panies gap-closing. Therefore the stable skyrmion texture of 
CMI derived from SHI does not carry nonzero quantum num- 
bers, which is consistent with the charge Chern number of the 
collinear magnetic ground state. 

To observe the inhomogeneity induced adiabatic responses, 
it is necessary to search systems in which both electron cor- 
relations and spin-orbit coupling are essential ingredients de- 
termining the material properties. Especially, to realize the 
TCMI, it is required to break lattice inversion symmetry. In 
this regard, the bilayer of perovskite-type transition-metal ox- 
ides grown along the [111] crystallographic axis is a promis- 
ing playground to realize the TCMI phase because the strength 
of the inversion symmetry breaking potential can be eas- 
ily manipulated in this system. The CMI with inhomoge- 
neous spin structures would be a fascinating venue to observe 
the interplay between the magnetism and electricity coupled 
through the nontrivial topological invariants. 

This work is supported by the Japan Society for the Pro- 
motion of Science (JSPS) through the "Funding Program for 
World-Leading Innovative R&D on Science and Technology 
(FIRST Program)". 

Appendix A: Adiabatic polarization current from Moyal 
product expansion 

In this Appendix, we explicitly demonstrate how the adi- 
abatic polarization currents of homogeneous/inhomogeneous 
systems can be derived using the gradient expansion (or 
Moyal product expansion) method. For the derivation, we 
have referred the discussions in Ref. l33l - [36l 

1. Polarization current in a system with translational 
invariance 

For a translationally invariant system, the total charge cur- 
rent along the a-direction {a = y) can be written as 

Jc^iej ^|^ti-[G-i(p,c^)9p„G(p,c^)], (Al) 

where the symbol tr denotes the summation over all discrete 
indices. Throughout this section, we keep the unit h because 
the gradient expansion is basically equivalent to the semi- 
classical expansion in powers of h. Let us first consider the 
adiabatic current induced by the smooth temporal variation of 
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the system. We assume that the system maintains the transla- 
tional invariance and the time dependence of the Hamiltonian 
appears due to the coupling of electrons to the time dependent 
the order parameter field N{t). Then the Green's function is, 
in general, a matrix possessing off-diagonal components in 
the time space such as G = G{p;t,t') = {t\G{p)\t'). The 
charge current can also be written as 



dttx{t\G-^{9)dp^G{9)\t). (A2) 



If the order parameter field N{t) changes slowly in time, 
the time dependence of the Green's function can be treated 
by using the Wigner transformation. In general, the Wigner 
transformation of a function Xi; ^2, X2) = F{xi\X2) is 
defined as 



where the center-of-mass coordinate R = ^^^"^ and the rel- 
ative coordinate r = x\ — with xx^i — (^i. 27X1.2). Also, 
■p — (w, p), p • r = p • X — bjt. If F is a convolution of two 
functions such as 



F(x,z)^ \ dyA{x,y)B{y,z). 



(A4) 



Then the Wigner transformations of F, A, B, which are repre- 
sented by F, A, B, respectively, satisfy the following relation, 

F{p,x) —A{p,x)B{p,x) 

-'-^[idpJ){dx.B) - {d^J){dpM 

-Y^idpA,A){d,^d,^B) + id,^d,J){dp^dp^B) 
- 2{dp^d^J){d^^dp^B)] + 0{h^). (A5) 



It shows that F can be written as a series of the space-time 
gradients of A and B in which the expansion parameter is 



given by h{^){-M-). In addition, the trace of F satisfies 



Tr[F] = / dxF{x, x) 



dxdp 



A{p,x)Bip,x), (A6) 



where d is the space-time dimension. 

Using the relation in Eq. (IA6I 1. the total charge current Ja 
(Eq. iA2\ ) is given by 



J Q, 



d^p f dtduj 



(27rft)2 J {2'Kh) 



tr[G(p,w,t)ap„if(p,w,t)], 
(A7) 



where K is the Wigner transformation of K = G ^ and G is 
the Wigner transformed Green's function that is defined as 

G(p,w,t) = y"d(ti-t2)e^"(*^-*^)G(p,ti,t2), (A8) 



where i = {ti +ti)l2. 

For the Hamiltonian li which is given by 



dp 



r7/-t(p,t)iJo(p,t)^(p,i), (A9) 



{2T^hY 

G and K — G^^ satisfy the following relations, 

ih-^ - Hoip, ti)] G(p, h,t2) = hS{ti - <2), (AlO) 

Oti J 



and 



J dt3K{p,ti,t3)G{p,t3,t2)^hS{ti~t2). (All) 



Therefore K is given by 

d 



Fip,R)= dre-^P'-FiR + r/2,R-r/2), (A3) K{p,ti,t2) 



ih—-Ho{p,h) 5{ti~t2) 
at I 



dn 

2TTh 



TT I ^1 + ^2 V 

0-i/o(p,^— ) 



e -iJ^(ti-t2)^ 
(A12) 

After the Wigner transformation of K, we obtain. 



K{p,u:,t) =uj- Haip,t) = g-\p,Lu,t), 



(A13) 



where Q is the semi-classical Green's function in which oj and 
t are treated as independent variables. 

The next step is to express the Wigner transformed Green's 
function G as a function of the semi-classical Green's func- 
tion Q. Let us first perform the Wigner transformation of 
Eq. ( lAl 11 1. From Eq. iA5i . we expect that the Wigner trans- 
formation of the left hand side of Eq. ( lAl li can be expanded 
in powers of h, which leads to following relation, 

00 

Q[K, G] = J2 fi^Qj G] = h, (A14) 
3=0 

in which 

eo[K,G] = KG, 

ei[x,G] 
e2[K,G] 



-[(ap.x)(a,.G)-(9,,if)(ap.G)], 

l[idpA,K){d.,^d^^G) + {d,,d^^K){dp,dp^G) 
o 

~2{dp^d^^K){d^^dp^G)], (A15) 



and so on. Now we make expansions of G and K in powers 
of h in such a way as. 



fc=o i=0 
Putting Eq. (IAI6I 1 into Eq. (IA14b . we obtain 

eo[i^o,fi^Go] = /i, 



(AI6) 



(A17) 
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and 



j+k+i>a 



(A18) 



where the summation is performed for j, k, I for given j + k + 
I > 0. Taking into account of the fact that K = Kq = G^^, 
Go, 1.2 can be obtained as 

Go = ^g, 

G2 = ^e[(5p.9p^.g-i)(5,.9,^.5) + {d..,d^,g-^)(dp,dp^g) 



2{dp^d,^g~^){d,^dp^g)] 



+ -r^g{d.,g-^){dp^ [g{dp^g-^){d.M) 



Ah? 
1 

1 



g{dp^g-^){d^^ [g{dp^g-^){d^M) 



(A19) 



Considering a smooth variation of the order parameter field 
N{t) in time, the leading order contribution to G is given by 



where v, A) = {uj, 9,px)- It is interesting to notice that the 
expression in the large parenthesis {} is nothing but the first 
Chern number, which is the topological origin for the quan- 
tized charge pumping through the cyclic variation of the order 
parameter field N{t). The electric charge induced by the one 
dimensional smooth structure can be described by the same 
topological invariant if t is replaced by the spatial coordinate 
X, as discussed in the recent work by Vayrynen et al.4^ 



2. Polarization current in a system with spatial inhomogeneity 

Now we consider the adiabatic polarization current in inho- 
mogeneous systems. In particular, we focus on the contribu- 
tion of the spatial inhomogeneity to the adiabatic current. For 
simplicity, we assume that the system maintains the transla- 
tional invariance along the x-direction and the spatial modu- 
lation occurs along the y-direction. Using the Wigner trans- 
formation, the total charge current of the system can be written 
as 



dy / dUi{t,y\G{px)dp^K{px)\t,y), 



2nh 

dpx f dydpy f dtduj 
2nh J 2TTh J 2TTh 

X tr[G{px;Py, y, w, t)dp^K{px;py, y; uj, t)] . 



(A23) 



G ^ng - —g[id^g-'){dtg) - idtg-^){d^g)]. (A20) 



Therefore, from Eq. ( IA7I ) and Eq. ( IA20b . the polarization cur- 
rent induced by the adiabatic temporal variation of the order 
parameter field N{t) is given by 



d^p f dtdu! 



{2TThY J (47rft) 

(dtg~^){d^g)}{dp^g~^)] 



iAQ{{dug-^){dtg) 



~eh 



2 j d^p f dtduj 

(27rfi)2 J {4:TTh) 



tv[{{d^g-')g{dtg-')g 



{dtg-')g{d^g-^)g}{dp^g-^)g], 

d^p f dtduj 



{2TThf J {l2Tlh) 

^e^,M{dt.g~^)g{dug~^)g{dxg-^)g], 

{All) 



where (/i, A) run over (oj, t,pa). In particular, considering 
a cyclic variation of N{t) = N{0{t)), which can be parame- 
terized using a polar angle 0(t) e [0, 27r], the adiabatic current 

along the a;-direction ji^"*, for example, can be written as 

^(1) _^ f dpy f f dpxdOdio 



2TTh 

X 



247r2 



e^.M{d^g-')g{d,g-^)g{dxg-')g]}, 

(A22) 



Smooth variation in the {t, y) space can be treated by using the 
gradient expansion method following the similar procedures 
taken in the previous subsection. 

We consider a class of Hamiltonians which can be written 

as 



H = 



^ dyij\px,y,t)Ho{px, -■^,y,t)'ip{px,y,t), 



{2TTh) 



(A24) 



in which Hq has the dependence on both y and . The in- 
verse Green's function K is given by 

K{px;yi,ti;y2,t2) 

d h d 1 

- Hq{px, -T—,yi,ti) S{ti - t2)S{yi - y2), 

oil I Oy^ i 



dfldp- 



- HQ{px,Py 



yi +2/2 ti+ t2 



{2nhf 
X giby(yi-2/2)-n(ti-t2)]_ 



After the Wigner transformation, K is given by 

K{px;py,y;uj,t) = w - Ha{px,Py,y,t) = g^^, 



(A25) 



(A26) 



where g is the semi-classical Green's function in which (a;, t) 
and {py, y) are treated as independent variables. 

To describe the inhomogeneity induced adiabatic current, 
we consider the gradient expansion of G in powers of h. Since 
we are interested in the leading order contribution of the spa- 
tial inhomogeneity to the adiabatic polarization current, only 
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the terms which are linear in the derivatives of t and y, re- 
spectively, are important. These are included in the G2 in 
Eq. ( IAI9I 1. The expression of the adiabatic current can be 
simplified by considering the following relations. At first, 
dp-dp.Q^^ — duidp^G^^ = from the general structure of 
the g\n Eq. (|A25ll.'' Similarly, du^dtG'^ = du^dyQ-^ = 0. 
Then the leading order contribution to the inhomogeneity in- 
duced polarization current can be summarized as 



t(2) _ j{2,A) , j{2,B) , j{2,C) , j{2,D) 



(A27) 



in which 



j(2,A) 



2TTh 



dydpy 



dtduj 



and 



j(2,S) 



4 J 2TTh J 2Trh J 2nh 

xti[g{dtdyg'')g{id^g-^)g{dp^g-')g{dp^g-') 
+ {dp^g-^)g{d^g-^)g{dp^g-^) 
+ idp^g-^)g{d^g-^)g{dp^g-^) 
+ idp^g-')g{dp^g-')g{d^g-') 
- id^g-')gidp^g-^)gidp^g-^) 
-{dp^g'')g{dp^g"')g{du.g'')}], 

(A28) 



ieh^ f dpx f dydpy f dtduj 



4 J 2nh J 2TTh J 2'Kh 

x\i[g{dtdp^g-^)g{idu.g-')gidyg'')gidp^g-') 
+ {dyg-^)g{d^g-^)gidp^g-^) 
+ {dp^g-^)g{d^g-^)g{dyg-^) 
+ idp^g-')g{dyg-')g{d^g-^) 
- {d^g-^)g{dp^g-')gidyg-^) 
^ {dyg-^)g{dp^g-^)g{d^g-^)}], 

(A29) 



and 



j(2,C) ^ 



hp f dpx f dydpy f dtduj 



4 J 2TTh J 2TTh J 2TTh 

{tr[g{dp^g-^)g{{dp^g-^)g{dx^g-^) 

+ {dx^g~^)g{dp^g-^)}g{dp^g-^)g{dx^g-')] 
+ tT[g{dx^g-')g{id.,^g-')g{dp^g-') 

+ {dp,g-')g{d,^g-')}g{dp^g-')g{dp^g-')] 
-tr[g{dx^g-')g{{dp^g-')gidp^g-') 

+ {dp^g-')g{dp^g-')}g{dp^g-')gidx^g-')] 
-tr[g{dp^g-^)g{{dx^g-^)g{dx^g-^) 

+ {d.,^g-^)g{dx^g-^)}g{dp^g'^)g{dp^g-')] 

(A30) 



and 



j: 



h^ f dpx f dydp. 



dtduj 



4 j 2TTh J 2TTh J 2Tih 

[ir[g{dp^g-')g{{dx^g-')g{dp^g-^) 

-{dp^g-^)g{dx^g-^)}g{dp^g-^)g{dx^g-^)] 
+ tx[g{dp^g-^)g{(dp^g'^)g{dx^g-^) 

-{dx^g-^)g{dp^g-^)}g{d.,^g-^)g{dp^g-^)\ 

(A31) 

j(2,A) , j(2,B). 



where {xi,pi) = {t,uj) or {y,py). Here Jx ' (Jx ' ) con- 



tains the second derivative of ^,i.e.,dtdyg ^ {dtdp g ^). 



Jx^''~'^ (Jx^'^^) consists of the terms which contains a sym- 
metric (antisymmetric) permutation of partial derivatives be- 
tween parenthesis in the middle. 

Lattice Dirac Hamiltonian: We apply the above result to the 
semi-classical lattice Dirac Hamiltonian, which has the fol- 
lowing structure. 



(A32) 



in which the five 4x4 traceless Gamma matrices sat- 
isfy {r^,r^} ~ 2(5^ with the identity matrix I. d — 
(di , (i2 , , (i4 , (is ) is a five-component vector that depends on 
the position x and momentum p coordinates. Using the iden- 
tity tr[rirjTfer;r„i] = ~4eijkim, it can be shown that 



(A33) 



(2 C) (2 D) 

Regarding Jx ' and Jx ' , it is crucial to observe the fol- 
lowing relation. 



j dLoty[g[d,,g-^)---g{d,,g-^)] 



dw- 
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' 1^1 1/2 "3 1^4 1^5 



tr[g{d,,g-') ■ ■ ■ g{d,,g-% 

(A34) 



where the five indices {1^1,1^2, 1^3, 1^4, 1^5) run over 
{u!,Px,Py,t,y). The antisymmetry of the trace under 
the permutation of the indices h'i leads to 



J 



(2,C) 



0, 



(A35) 



and 



Ji 



2,D) 



ieh^ f dpx f dydp. 



dtduj 



60 J 2TTh J 2TTh J 2TTh 

X .,.3.4.5 tr[^(a.,^-i) • • • gid,,g-')] 



(A36) 



For the description of the charge pumping in topological 
chiral magnets, we can assume that the spatial inhomogeneity 
of the system appears only through the coupling to the inho- 
j I mogeneous order parameter field N{y, t) and the charge cur- 
rent is induced by the adiabatic variation of N{y, t), which 
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satisfies N"^ = 1, i.e., N € S"^ . We introduce a spherical 
coordinate (6', for the space-time (y, i) dependent unit vec- 
tor field N{y,t) = N{e{y,t),(j){y,t)). Then J^^'"^ can be 
rewritten as 

=e4'hL. (A37) 

where the indices (/ii, ^2, MSi /^47 /^s) run over 
{uj,Px,Py,0,(t)). Here indicates the 2-torus given by 

two dimensional Brillouin zone. It is interesting to note 

(21 

that Cf^jig^^ is nothing but the second Chern number, which 



is quantized when the spherical solid angle $7(6', 0) is fully 
covered by N{y,t). Therefore the inhomogeneity induced 

adiabatic current J-i can be summarized in the following 
suggestive way, 

7(2) = 7(2) , 7(2) .AC!0>) 

'^x "^x, topological ' "^x, perturb' i,J-vjoy 

(2) (2 D) 

m which topological = Jx IS the topological contri- 
bution to the adiabatic current whose topological nature is 
endowed with the second Chern number All the other 
non-topological (perturbative) contributions are included in 

.(2) _ j{2,A) j(2,B) j{2,C) 

"^x.peiturb — + Jx + Jx 
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